Graphene field effect transistors with ferroelectric gating 
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Recent experiments on ferroelectric gating have introduced a novel functionality, i.e. non- 
volatility, in graphene field effect transistors. A comprehensive understanding in the non-linear, 
hysteretic ferroelectric gating and an effective way to control it are still absent. In this letter, we 
quantitatively characterize the hysteretic ferroelectric gating using the reference of an independent 
background doping (ubg) provided by normal dielectric gating. More importantly, we prove that 
«bg can be used to control the ferroelectric gating by unidirectionally shifting the hysteretic ferro- 
electric doping in graphene. Utilizing this electrostatic effect, we demonstrate symmetrical bit writ- 
ing in graphene-ferroelectric FETs with resistance change over 500% and reproducible no-volatile 
switching over 10 cycles. 
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The electric field effect, which continuously tunes the 
Fermi level (Ep) in the conical energy band structure of 
graphene, plays a critical role in studying the extraordi- 
nary electronic properties of graphene [l| ■ Using conven- 
tional dielectrics such as Si02 and more recently Hf02, 
PMMA and AI2O3 with linear dielectric response to elec- 
tric field, many fascinating physics have been discovered. 
Among these celebrated phenomena are the anomalous 
quantum Hall effect 0, Klein tunneling and gate- 
tunable bandgap in bilayer graphene 4]. Despite such 
tremendous progress, there is a keen interest in the sci- 
ence community to utilize new dielectrics and substrates 
for exploring new graphene physics and functionalities 
P-Q. Among promising candidates, ferroelectrics are 
unique both in ultra high dielectric constants (k) up to 
a few thousands and non- linear, hysteretic dielectric re- 
sponse to electric field. The ultra high n makes ferro- 
electrics promising substrates for studying charge scat- 
tering mechanism in graphene which could be a 
crucial step in realizing ultra-high mobility 14j in non- 
suspended graphene. Equally important, the ultra-high k 
may allow ultra high doping in graphene with charge den- 
sities (> 10 14 cm -2 ) exceeding electrolyte doping [15j and 
with gate tunability at cryogenic temperatures. Based on 
the hysteretic ferroelectric gating, a novel functionality of 
non-volatile graphene-ferroelectric field effect transistors 
(GFeFETs) has been demonstrated 

However, the fundamental understanding of ferroelec- 
tric gating is still elusive. In contrast to the linear doping 
vs normal dielectric gating relation, n = aV g \^, ferro- 
electric gating introduces a pronounced hysteresis in the 
charge doping. In particular, polymer ferroelectric gat- 
ing introduces strong electron-hole puddles in graphene 
even far away from the Dirac point. Therefore, Hall mea- 
surements alone may be misleading in determining the 
induced charge doping. Thus, a quantitative modeling 
will not only improve the understanding of ferroelectric 
gating but also help in optimizing the performance of 
GFeFETs. Ferroelectric gating is also characterized by 
two symmetrical remnant polarizations, i.e. P-j- = P r and 



P^ = — P r for upwards and downwards dipole configu- 
rations, respectively. Consequently, P-j- and Pj. induce 
two identical zero-field resistance states in graphene. Al- 
though two distinct resistance states can be created by 
polarizing (Ro(P r )) and depolari zing (Ri(P ~ 0)) the fer- 
roelectric thin film alternately 16|, the depolarization 
state is not in thermodynamic equilibrium and less sta- 
ble than the polarization state. To solve this problem, 
we need an effective way in controlling the hysteretic fer- 
roelectric doping. Last but not least, GFeFETs in our 
earlier work [16| is characterized by low charge carrier 
mobility of few hundred cm 2 V _1 s _1 . Such low mobil- 
ity prevents the determination of the intrinsic physical 
properties and limitations of GFeFETs. 

In this letter, we present a quantitative understand- 
ing of high quality graphene devices under ferroelectric 
gating. For this purpose, we introduce an independent 
reference doping (tibg) by the Si02 back gating. We 
show that the evolution of the device resistance hystere- 
sis from symmetrical double peak to asymmetrical single 
peak structures can be consistently simulated by the elec- 
tric displacement continuity equation using the reference 
of the Si02 gating. We also show that by controlling the 
polarity and magnitude of ubGi the hysteretic ferroelec- 
tric doping in graphene can be shifted unidirectionally. 
In analogy to exchange biased spin valves [l7|, this ef- 
fect provides a reference point for maximizing the resis- 
tance change at zero electric field and enables symmet- 
rical bit writing in GFeFETs. We demonstrate highly 
reproducible non-volatile switching over 10 5 cycles and 
AR/R exceeding 500% in GFeFETs. 

The GFeFET sample geometry is shown in Fig. QJi. 
Detailed sample fabrication procedures have been dis- 
cussed in Ref. [16|. For the samples used in this study, 
the ferroelectric thin film of poly(vinylidene fluoride- 
trifluoroethylene 72:28) (PVDF) is ~ 0.5 /im thick. The 
GFeFETs were electrically characterized at room tem- 
perature in vacuum using four-contact lock-in technique. 

Before polarizing the ferroelectric, we first measured 
the Hall mobility and the resistance vs Si02 gate voltage 
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FIG. 1. (a) Sample geometry of GFeFETs. (b) R vs Vbg of 
one sample after PVDF coating. Red open square and black 
solid line are the experimental and fitting results respectively. 
Inset: Atomic force microscopy of the sample after PVDF. 
Color scale: to 164 nm. 



characteristics (R vs Vbg) to determine the sample qual- 
ity. Most samples retain their high mobility after PVDF 
spin-coating and annealing, as shown in Fig. [T|b for a typ- 
ical sample with Hall mobility of 4,600 cm 2 V _1 s _1 (18| . 
Quantitatively, the ambipolar R vs Vbq characteristics 
can be fitted very well by the model [19(, 
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using the Hall mobility fiHaU- For the sample shown in 
Fig. QTj, the fitting yields a residual carrier concentration 
?i rcs = 2.77 x 10 n cm- 2 . 

Compared to the Si02 gating, one fundamental differ- 
ence introduced by ferroelectric gating is pronounced hys- 
teresis in the resistance vs ferroelectric gate voltage char- 
acteristics (R vs Vtg)- Though such hysteretic R vs Utg 
can be qualitatively explained by the electric displace- 
ment continuity equation at the ferroelectric/graphene 
interface [l6j |. a quantitative understanding of ferroelec- 
tric gating is still missing. Here, we introduce an in- 
dependent «bg using the SiC^/Si back gate. This pro- 
vides a well defined, constant reference for determining 
the doping induced by PVDF gating. To study the effect 
of tibg on the ferroelectric gating of GFeFETs, it is also 
important to limit the polarization magnitude in PVDF, 
since the effect of ferroelectric gating is nearly 10 times 
stronger than the Si02 gating [16(. Thus, we first intro- 
duced very small |P r | in PVDF by limiting the maximum 
top gate voltage (UrGmax) to ±5 V. Such low Vtg only 
slightly polarizes PVDF, allowing tibg to match or even 
exceed the |P r | induced doping in graphene. 

In Fig. (2J}, we show the resistance of the GFeFET as 
a function of both Vtg and Vbg- With Vbg ~ 6 V, the 
R vs Vtg curve shows two symmetrical resistance peaks 
and nearly negligible AR/R (Fig. [2p(iii)). By gradu- 
ally tuning hbg with Vbg, the two resistance peaks be- 
come more asymmetrical and shift leftward (rightward) 
for nee < (tt-bg > 0). The shift in peak positions 
leads to an increase in AR/R, which has a maximum at 



V BG » -6 V (Fig. [2p(ii)) and Vbg « 18 V (Fig. [2p(iv)), 
respectively. Crossing these two points, AR/R decreases 
as |tibg| keeps on increasing. At large enough 7Jbg, the 
double peak structure eventually disappears in the R vs 
Vtg hysteresis (Fig. [2b (i) and Fig. db(v)). 

The evolution of the resistance peaks and the change 
in AR/R can both be explained by two independent but 
competing doping processes in graphene by polarized fer- 
roelectric dipoles and Vbg , respectively. For such a dual- 
gated system, the interfacial electric displacement conti- 
nuity equation is expressed by 



/ 8P(V T G) + n* = n(V T G,VBG)e, 



(2) 



where /3P(Vtg) represents the hysteretic dipole doping 
by the ferroelectric gating [2(J, and n* = n cm + tibg 
is the reference doping induced by the dielectric envi- 
ronment and Vbg respectively. For n* s» 0, the doping 
in graphene is dominated by the ferroelectric gating by 
n(Vr G ,n* « 0) = -/3P(U TG )/e. Using Eq. [TJ it is now 
straightforward to see that n(VrG,n* rj 0) will produce 
a R vs Vtg hysteresis with two symmetrical resistance 
peaks, centering on the two coercive- field points where 
P(Vtg) crossing zero. Experimentally, this is the R vs 
Vtg curve in Fig. [2]b (iii) with Vbg = 6 V, in which two 
resistance peaks are centered at Vtg = ±2.2 V respec- 
tively. By converting each R in Fig. [2jD(iii) into doping 
using Eq. [TJ we directly determined the doping curve 
n(VrG,«* ~ 0). The result is shown in Fig. [2fc (red 
curve). As expected, this doping curve is hysteretic and 
characterized by two zero-field doping levels with equal 
magnitude, i.e. \n\ \ = |no| = /3P r /e. 

After acquiring n(VrG, n * ~ 0), we can deduce indi- 
vidual R(Vtg>"*) curves for non-zero n* by substituting 
n( Vtg, U BG ) = -pP(V TG )/e + n em + aV BG into Eq. [TJ 
Here a = 7.2 x 10 10 cm" 2 V _1 is the doping coefficient 
of 300 nm Si02, and fitting parameter [2l| . 

By tuning n env and matching the resistance peaks of the 
simulation to the experimental, we simulated each exper- 
imental R(VtG) Ubg) curve in Fig. 2b. As shown by blue 
dotted lines, the simulation reproduces the evolution of 
the experimental results very well. Two resulting dop- 
ing hysteresis for the resistance curves in Fig. 0)(i) and 
Fig. [2Jd (ii) are further compared with n(VrGi n * ~ 0) m 
Fig. [2t. From the comparison, we can see that AR/R 
approaches the maxima as one zero-field doping level sits 
near the Dirac point when |n*| ~ /3P r /e (blue hysteresis 
loop). Further increase in tibg moves both n\ and no 
away from the Dirac point, and AR/R decreases (black 
hysteresis loop). 

Thus, we have shown that using a background dop- 
ing introduced by normal dielectric gating as a refer- 
ence, the hysteretic behavior of R vs ferroelectric gating 
in GFeFETs can be quantitatively determined by solv- 
ing the electric displacement continuity equation. For 
memory applications, AR/R is of great importance. Fol- 
lowing the above discussions, the two zero- field resis- 
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FIG. 2. (a) R vs Vtg and Vbg of the GFeFET with very small |P r |. (b) Extracted single traces of R vs Vtg with different 
Vbg- The blue dotted lines are simulated results, (c) ubg tunable doping hysteresis in GFeFETs. Only three doping hysteresis 
loops, corresponding to experimental curves i, ii and iii in Fig. [2}j, are shown for clarity. See main texts for discussions. 
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FIG. 3. (a) AR/R as a function of Vbg with different VTGmax. 
Two maxima are observable with VTGmax = 5 V (black open 
circles). The red solid line shows the simulation with /JP r = 
4.2 x 10 11 cm" 2 . For V T Gmax = 30 V, the maximum AR/R 
is increased to 500% (Vbg = 32 V). (b) R (Vtg, Vbg) of the 
GFeFET with higher /3P r (~2x 10 12 cm" 2 ). Double peak 
structures dominate over the whole Vbg range. 



tance states are Ri = , and Rn = 

V7e^^/n? oa + (/3P r /c-n*) 2 

; L respectively. Thus, the best strat- 

egy to utilize the field-dependent resistance is to fully 
polarize the ferroelectric and introduce a matching ubGi 
as demonstrated in Fig. With VrGmax = 5 V, two 
maxima of ^250% are present in AR/R vs Vbg, which 
can be also simulated very well by Eq. Q] and [2] with 
/3P r /e = 4.2 x 10 n cm" 2 . By increasing VTGmax to 30 
V, the maximum AR/R is increased to 500%. The fast 
increase in P r not only increases the maximum AR/R, 
but also increases the separation between the two AR/R 



maxima, resulting in one maximum being outside of Vbg 
measurement range. For this VrGmax , R. vs Vtg shows a 
dominant double peak structure over the full Vbg range 
(Fig. [3Jd). However, we can still see the tendency of 
a transition from double peak structure to single peak 
structure as Vbg exceeding 40 V. 

Such rtBG shifted hysteretic doping in graphene is a 
ferroelectric analogy to the ferromagnetic exchange bias 
17j . Utilizing this electrostatic effect, the bit writing in 
GFeFETs can be much simplified by switching the ferro- 
electric polarization between P r and -P r , using symmet- 
rical voltage sweeps. With tibg ~ ~ /?P r /e, to write the 
high resistance "1", a negative writing voltage (-^writing) 
is applied to the ferroelectric, setting the dipole polariza- 
tion to -P r independent of the initial states in the unit 
cell (Fig. andSh). In contrast, a positive Vwdting with 
the same magnitude sets the GFeFET into low resistance 
"0" (Fig. 2J: and[4ji). Compared to the asymmetrical bit 
writing by polarizing and depolarizing the ferroelectric 
alternately [l6(, such symmetrical writing in GFeFETs 
not only provides simplicity but also takes full advantage 
of the fast switching speed of ferroelectric. For lead zir- 
conate titanate (PZT) based materials, this can be as fast 
as 280 ps . Another potential application of this elec- 
trostatic effect could be multi-bit-per-cell data storage in 
GFeFETs utilizing the ti B g tunable AR/R. 

We have also tested the reproducibility of our GFe- 
FETs working with j3Y I rj \n*\e. During the fatigue test, 
a triangular wave of lk Hz was applied to the PVDF thin 
film. Every 12 (24) seconds, the triangular wave was in- 
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FIG. 4. Symmetrical bit writing in GFeFETs with tlbg ~ 
—/3P T /e. (a) and (b) Writing "1" using — Uwriting- (c) and 

(d) Writing "0" using Venting- Dashed and solid arrows in- 
dicate the forward and backward voltage sweep directions re- 
spectively. The writing is independent on the initial states. 

(e) and (f) Fatigue test of one GFeFET with symmetrical bit 
writing, showing non- volatile switching cycles exceeding 100k. 



terrupted and one R vs Vtg curve was recorded. The 
corresponding AR/R as a function of switching cycles 
and the raw data of individual R vs Vtg curve are sum- 
marized in Fig. and |3f respectively. The fatigue test 
clearly demonstrates reproducible non-volatile switching 
exceeding 100k cycles in the GFeFET. Ultimately the 
life span of PVDF-based GFeFETs is 10 7 0. Thus, 
PVDF-GFeFETs could provide a cost-effective solution 
for flexible non- volatile data storage with sub-/xs switch- 
ing speed. On the other hand, inorganic ferroelectric 
(such as PZT) should be used if fast writing speed (< 
ns) and ultra-high endurance (10 10 ) are required. 

In conclusion, we have demonstrated a quantitative 
way in determining and controlling the hysteretic ferro- 
electric gating in GFeFETs. Using an i ndep endent linear 
dielectric gating (tibg) as a reference [2J], ferroelectric 
gating can be quantitatively determined by the electric 
displacement continuity equation. The reference gating 
can also be used to control ferroelectric gating by intro- 
ducing a unidirectional shift in the hysteretic ferroelec- 



tric doping in GFeFETs. One specific application of this 
electrostatic "bias" effect is symmetrical bit writing in 
GFeFETs directly utilizing P r and — P r with speed and 
simplicity. The ferroelectric gating phenomena and re- 
lated modeling and controlling methods presented in this 
study will be important in understanding future charge 
transport studies on ferroelectric gated graphene elec- 
tronic devices. 
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